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Starting from the simplified analytic model of electronic spectrum of iron - pnictogen (chalcogen) high - 
temperature superconductors close to the Fermi level, we discuss the influence of antiferromagneting (AFM) 
scattering both for stoichiometric case and the region of possible short - range order AFM fluctuations in 
doped compounds. Qualitative picture of the evolution of electronic spectrum and Fermi surfaces (FS) for 
different dopings is presented, with the aim of comparison with existing and future ARPES experiments. 
Both electron and hole dopings are considered and possible pseudogap behavior connected with partial FS 
"destruction" is demonstrated, explaining some recent experiments. 



'PACS: 71.10.Hf, 7f.l8.+y, 74.25.Jb, 74.70.-b, 74.72.-h 

Recent discovery of the new class of iron based high- 
|temperature superconductors [1] stimulated intensive 
of experimental and theoretical efforts to understand 
lits properties (see for the review Refs. [H I2)- De- 
spite aheady the immense progress in understanding 
of these systems, the nature of superconducting pair- 
ing and anomalies in the normal state are still under 
debate. 

Clarification of the structure of electronic spectrum 
of new superconductors is crucial for explanation of 
.their physical properties. Accordingly, since the first 
days, different groups have started the detailed band 
structure calculations for all classes of these com- 
pounds, based primarily on different realizations of gen- 
eral LDA approach. These calculations were primarily 
performed for paramagnetic tetragonal FeAs 1111 sys- 
tems [lISlinilT], for 122 [iHIig, for 111 [inillTlIIl] and 
|a-FeSe [T3], followed by many similar works by other 
authors. In fact, all these calculations demonstrated al- 
,most universal LDA band structure in relatively narrow 
'energy interval (iO.ley) around the Fermi level, which 
lis of relevance to superconductivity [2] . 

In this energy interval the electronic spectrum can 
be modelled analytically as follows. Three "hole-like" 
^branches of the spectrum crossing the Fermi level near 
the F point in the Brillouin zone (cf. FigH^) can be 
taken isotropic and modelled by quadratic dispersion: 






(1) 



where rrii, Si (i — 1,2,5) can be easily determined from 
LDA calculations (e.g. for 122 system from the results 
ofRef. 0). 

Two "electron-like" branches of the spectrum cross- 
ing the Fermi level near M(7r, tt) point of the reduced 
Brillouin zone are anisotropic and produce two el- 
liptic isoenergetic crossections at the Fermi level (cf. 
FigUlD), one of which is extened in the direction MF, 
with the second one extended in the orthogonal direc- 
tion. Let us count the momentum from the M point 
(i.e. replace p — Q — > p) and take one momentum 
p axis along MF direction and other orthogonal to it 
(FiglTJa). The relevant momentum projections pi and 
P2 are connected with the usual x, y projections as 
Pi = ^-^h§^, P2 = ^^~^'' . Consider one of the eUipses, 
e.g. those extended along the direction orthogonal to 
MF direction. Electron dispersion along MF can be 

2 

modelled by quadratic law £pi(p) — -^^ — £4. Dis- 
persion along the orthogonal (to MF) direction is de- 
termined by higher (in energy) branch of the spectrum, 
originating from the hybridization of two "bare" disper- 
sions (cf. Fig. [T^), which we also assume quadratic. 

Then, neglecting the small hybridization gap, we ob- 

2 2 

tain ep2ip) = Max{-^ - e^,^ - £5) Parameters 

TO4, £4, TO5, £5, Too can be taken from LDA data. Thus, 

for anisotropic "electron-like" spectrum we can use the 

following model: 



£4(p) 



\cb)epi{p) + stn\^)ep2{p) (2) 
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where p^ = Pi + P2 1 ^^id (p is the polar angle with re- 
spect to pi axis. This model guarantees correct energy 
crossections in direction MF and orthogonal it, as well 
as isotropy of the spectrum in case of £pi (p) — £p2 (p) ■ 
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Fig.l. Qualitative picture from of the band structure in 
Mr direction in the reduced Brillouin zone (a) and the 
relevant Fermi surfaces (b). 

Energy dispersion for the second "electron-like" band 
e5(p) is also given by Eq. ^ virith the obvious substitu- 
tion -H> f + '/'• Finally, we describe the "electron-like" 
bands in our model as: 



P2 



£4(p) = <^ ^p; 2"^" 
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-£4 



for p2 = pl+pj < p2 



P2 Pi. 



£4-^65 for p^ > pt, 



(3) 

where p^ — 2{e^ ~ ^i)/{^ ^ ^) i^ the square of the 
momentum at the crossing of two "bare" hybridizing 
bands. 

The qualitative picture of electronic spectrum and 
Fermi surfaces is shown in Fig. [TJ Essentially this 
kind of electronic spectra and Fermi surfaces in new 
superconductors were qualitatively confirmed by angle 
resolved photoemission spectroscopy (ARPES), starting 
with the early works [H [IB [H [Q III [H [23 IH] , fol- 
lowed by many further studies by the same and other au- 
thors. Most of these experiments were performed on sin- 
gle crystals of 122 systems, while for other compounds 
good quality single crystals are up to now just unavail- 
able. Though in qualitative agreement with the results 
of LDA calculations, these experiments show rather dif- 
ferent results concerning finer details, such as the pre- 
cise number of "hole-like" FS cylinders around the F 
point, as well as the topology of "electron-like" cylin- 
ders around the M - point. 

In general LDA calculations underestimate the role 
of electronic correlations. ARPES experiments show 
that these systems apparently belong to the class of 
intermediately correlated systems, with correlation in- 
duced band narrowing by the factor of two [16| . This 
is confirmed by some of LDA-f DMFT calculations [35] , 
though theoretical situation here remains rather con- 
troversial. In the following we take correlations into 
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Fig. 2. Recurrence "Dyson equation" for the Green's 
function. 



account by simple rescaling of the energy by the factor 
of two as compared with LDA [TB] . 

Undoped FeAs compounds are antiferromagnetically 
ordered with AFM vector Q — (0, tt) in extended Bril- 
louin zone, corresonding to Q = (tt, tt) in the reduced 
zone [21 13] . Electron or hole doping suppresses AFM or- 
dering and induces superconductivity, similar to the well 
known situation in cuprates. Recent neutron scattering 
experiments [531 HI] clearly show that in the substan- 
tial part of the phase diagram of FeAs systems in normal 
paramagnetic state rather strong fluctuations of AFM 
short-range order persist, as predicted e.g. by the model 
of "nearly antiferromagnetic Fermi liquid" [25] [26] [22] ■ 
These fluctuations can, in principle, induce the pseu- 
dogap behavior in electronic spectrum, similar to that 
observed in cuprates [28j . 

Effective interaction of electrons with AFM spin fluc- 
tuations is determined in this model by dynamic spin 
susceptibility characterized by the maximum at scat- 
tering vectors close to AFM vector Q = (tt, tt), which 
we assume here to be the same for electron from differ- 
ent bands and for interband scattering. Limiting our- 
selves to high enough temperatures we can neglect the 
dynamics of AFM fluctuations and consider them Gaus- 
sian [3S]. The Green's function for electrons moving in 
the "quenched" Gaussian random field of these fluctu- 
ations can be represented by recurrence "Dyson equa- 
tion" shown in Fig[51 which is the direct multiple bands 
generalization of the summation procedure, proposed 
and actively used in Refs. [IH] 133 Ell 132] > taking into 
account all Feynman diagrams for electron scattering in 
such random field. 

Analytically, this "Dyson equation" can be written 
as: 

Gl = GZ,5,, + G-,A'<n + 1) ^ G^+/ ^ G\) (4) 

kfn I 

where i,j represent band indices, A characterizes the 
AFM pseudogap width (of the order of AFM band split- 
ting). 



G??.(i?p) = 



1 



(5) 



E — e"(p) + inv^K. 

K = £^^^ is an inverse correlation length of AFM short- 
range order fluctuations, ef (p) = ei(p + Q) and w" = 
Iff (p -I- Q)|-H|wf (p + Q)| for odd n, while ef (p) = e^(p) 
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and w" = l^f (p)l + bi'lp)! for even n. Velocity projec- 
tions wf (p) and wf (p) are determined by the momen- 
tum derivatives of electronic dispersion in the i-th band 
£i(p). Combinatorial factor s{n) for the case of Heisen- 
berg AFM fluctuations (spin-fermion model of Ref. [31]) 
is given by: 



s{n) 



2d^ for odd n 
^ for even n. 



(6) 



The physical Green's function corresponds to n = 0. 
Then, after some simple manipulations we may show 
that 

where the physical self-energy 

E(^p) = E"=i(^p) 



(8) 



is determined from the recurrence procedure (continued 
fraction representation) : 



E"(i?p) = 



A^s{n) 



(Go"(£;p))-i-S"+i(i^p) 



(9) 



where Gq{Ep) ~ J2j Gqj{Ep). As a byproduct of these 
general equations we can easily analyze the electronic 
spectrum in the case of AFM long-range order, trun- 
cating the continuous fraction in Eq. (j9]) at n = 1 and 
taking the limit of k — > 0. The spectral density and 
density of states are obviously given by: 



A{Ep) 



-ImSpG^iiEp)- N{E) 

IT 



P 



A{E,p) 



(10) 

We performed calculations for a variety of parame- 
ters of the model, using for the spectrum LDA data for 
122 from Ref [5], scaled by factor of two to account for 
correlations. Below we present results for A = 50 meV, 
which is in rough agreement with the estimates of AFM 
band splitting from ARPES data [331 131] and neutron 
scattering [35 (varying in the interval 50-100 meV), cor- 
relation length of AFM fluctuations ^ = 10a (a - lattice 
spacing) , also in rough agrrement with netron scattering 
data [23l [24] . In the following, all momenta are given 
in units of inverse lattice spacing, energies in eV. To 
make the results comparable with ARPES experiments 
we have also introduced effective widening to simulate 
finite energy resolution of ARPES replacing E -^ E + i'y 
and taking 7=10 meV (corresponding to best ARPES 
resolution). 

In Fig. [3] we show "ARPES" energy bands of 122 
system, revealed by the maps of spectral density, along 



main symmetry directions, starting from the case of nor- 
mal (paramagnetic) LDA bands, via AFM long-range 
ordered state, to "pseudogapped" state, characterized 
by electrons scattered by short-range ordered AFM fluc- 
tuations - AFM band splittings transforming to pseu- 
dogaps due to AFM short-range order. 

In Fig. [H we show spectral density maps at the 
Fermi level for different dopings - from slightly elec- 
tron doped, via undoped, to hole underdoped and opti- 
mally hole doped case. These maps essentially produce 
"ARPES" Fermi surfaces of 122 system at different dop- 
ings. In fact, the system always remains metallic in a 
sense that at every doping we observe "open" Fermi sur- 
faces, though we also can see rather complicated series of 
Fermi surface transformations, with some cylinders be- 
ing almost "destructed" (damped) either by AFM long- 
range order, or by short-range order AFM fluctuations. 
Of these maps, we identify the last one in the third 
row (4c) as corresponding more or less to optimally hole 
doped case in satisfactory agreement with ARPES data 
e.g. from Refs. [HI [T51 [501 [35] , while the third one in 
the same row (3c) apparently well corresponds to the 
hole underdoped case studied in Ref. [36J, demonstrat- 
ing the inner hole cylinder rather damped by pseudo- 
gap fluctuations with characteristic wave vector of the 
order of AFM vector Q. Significant pseudogap forms in 
the (partial) density of states on precisely this cylinder, 
in agreement with Ref. |M|. In general, the available 
ARPES data suffer from rather bad resolution, so that 
pseudogap fluctuations can significantly complicate ob- 
servation of all FS cylinders and much work is needed 
to reveal possible complicated picture of FS transfor- 
mations, illustrated in Fig. [U It should be taken into 
account that pictures shown in the second row (b) of 
Fig. [Hare sensible only within the part of the phase di- 
agram with AFM long-range order, while the third row 
(c) applies to paramagnetic region, where superconduc- 
tivity appears at lower temperatures. 

Our calculations show, that the pseudogap forms 
only in (partial) densities of states, corresponding to 
those cylinders strongly affected by short-range AFM 
fluctuations, and this is not, in general, "pinned" at the 
Fermi level. Pseudogap in the total density of states is 
always rather weak, and the problem remains, whether 
it is sufficient to explain claims for the pseudogap be- 
havior observed in some NMR experiments [3J [^ . 
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Fig.3. Energy bands. Upper panel (a) - "bare" (scaled LDA) bands in paramagnetic state in the absense of AFM 
fluctuations. Panel (b) - AFM long-range ordered state with A =0.05 eV. Panel (c) - bands in the pseudogap state 
induced by AFM short-range order fluctuations with ^ = IQa and A =0.05 eV. All bands are shown with finite 
"experimental" resolution 7 =0.0f eV. Dotted lines show Fermi levels for different dopings used in our calculations 
of Fermi surfaces below. 
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Fig.4. "ARPES" Fermi surfaces at different doping levels shown by dotted lines in Fig[3] Column 1 - electron 
doping with Ef =0.02 eV, 2 - undoped system with Ef =0, 3 - hole doping with Ef =-0.035 eV (hole underdoped 
system), 4 - optimal hole doping with Ef =-0.085 eV. Upper panel (a) - "bare" FS in paramagnetic state in the 
absense of AFM fluctuations. Panel (b) - AFM phase with A =0.05 eV. Panel (c) - pseudogap state with ^ =10a 
and A =0.05 eV. 



